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Abstract 

A formalism based on a relativistic plane wave impulse approximation is developed to investigate 
the strange-quark content (g A ) of the axial- vector form factor of the nucleon via neutrino- nucleus 
scattering. Nuclear structure effects are incorporated via an accurately calibrated relativistic mean- 
field model. The ratio of neutral- to charged-current cross sections is used to examine the sensitivity 
of this observable to g s A . For values of the incident neutrino energy in the range proposed by the 
FINeSSE collaboration and by adopting a value of g s A = —0.19, a 30% enhancement in the ratio is 
observed relative to the g% = result. 
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I. INTRODUCTION 



In a desperate effort to salvage the conservation of energy and angular momentum, Wolf- 
gang Pauli postulated the existence of the neutrino: a neutral and (almost) massless particle 
that he feared will never be detected. Seventy five years later and with three neutrino fami- 
lies firmly established, neutrino physics in the 21st century has taken center stage in fields as 
diverse as cosmology, astro, nuclear, and particle physics. With the commission of new neu- 
trino observatories and facilities for the study of solar, atmospheric, and reactor neutrinos, 
conclusive evidence now exists in favor of neutrino oscillations. Although neutrino-oscillation 
experiments have now evolved from the discovery into the precision phase, important ques- 
tions remain unanswered pj. Depending on these answers, a radical modification to the 
standard model of particle physics may be required. At the very least, neutrino oscilla- 
tions already demand a mild extension to the standard model: in the standard model the 
individual lepton numbers must be conserved. 

The discovery of neutrino oscillations established two incontrovertible facts: a) that neu- 
trinos have mass and that these masses can not all be equal and b) that the three known 
neutrinos z/ e , u^, and v T are linear combinations of three neutrino mass eigenstates (com- 
monly referred as u%, i>2, and z/3). Note that due to the very small neutrino masses, os- 
cillation experiments are sensitive to only the squared mass difference (Am 2 ) of the mass 
eigenstates. Further, the unitary matrix linking the flavor to the mass eigenstates has a 
well-known counterpart in the quark sector: the CKM matrix. While precision experiments 
have started to determine squared mass differences and mixing angles, a strong, vibrant, 
and interdisciplinary program is now being established to tackle the myriad of remaining 
open questions pj. Among these are: What is the mass hierarchy? Is the neutrino a Dirac 
or Majorana particle? Does the neutrino matrix contain a CP-violating phase that may 
explain our matter-dominated Universe? Is there a need for additional "sterile" neutrinos? 

The apparent need for an additional sterile neutrino stems from two older experiments. 
The first one in 1989 measured the total and partial (into hadrons and charged leptons) 
widths of the Z° boson at the large electron-positron (LEP) collider at CERN and extracted 
the number of neutrinos flavors to be N v = 3.00 ± 0.08; note that a more recent analysis 
reports N v = 2.984 ± 0.008 0. The second experiment was a 1995 neutrino-oscillation 
experiment at the Liquid Scintillator Neutrino Detector (LSND) at the Los Alamos National 
Laboratory [3J. This experiment reported evidence of — > v e oscillation, but with a squared 
mass difference that is too large — and thus inconsistent — with the two independent values 
extracted from solar, atmospheric, and reactor experiments. Simply put, for three neutrino 
flavors the algebraic relation Am^ + Am^Am^ must be satisfied, yet Am^ + Am^^ 
^ m LSND- The most favorable scenario that accommodates three independent Am 2 values 
is the addition of a sterile neutrino. Yet other possibilities exist: the LSND analysis may be 
incorrect. It is the primary goal of the Booster Neutrino Experiment (B00NE and its first 
phase MiniBooNE) at Fermilab to confirm the LSND result. 

While MiniBooNE's primary goal is to confirm the LSND result, this unique facility is also 
ideal for the study of supernova neutrinos, neutrino-nucleus scattering, and hadronic struc- 
ture. An ambitious experimental program — the Fine-grained Intense Neutrino Scattering 
Scintillator Experiment (FINeSSE) — aims to measure the strange-quark contribution to 
the spin of the nucleon |H H E| • FINeSSE is part of a larger program started in the late 80's 
that attempts to answer a fundamental nucleon-structure question: how do the non- valence 
( "sea" ) quarks — particularly the strange quarks — contribute to the observed properties 
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of the nucleoli? To date, the role of strange quarks in the nucleon remains a contentious 
issue and one that remains a subject of intense activity all over the world. In an attempt 
to find a satisfactory answer to this fundamental question, a number of reactions have been 
proposed. These include: (i) deep inelastic scattering of neutrinos on protons 0, H|, (ii) 
deep inelastic scattering of polarized charged leptons [9f, (iii) pseudoscalar meson scattering 
on a proton and (iv) parity- violating electron scattering 0, Part of the contro- 



versy arises because these reactions do not all reach similar conclusions. For example, both 
reactions (i) and (ii) suggest a non-zero strangeness contribution, in contrast to reaction 
(iv) that indicates a strange quark contribution to the charge and magnetic moment con- 
sistent with zero. Parity-violating electron scattering in particular has received extensive 
experimental attention as in the SAMPLE Collaboration at the MIT-Bates accelerator , 
the HAPPEX Collaboration at the Jefferson Laboratory and the A4 Collaboration at 
the MAMI facility in Mainz 0|. However, theoretical investigations have shown that large 
radiative corrections and nuclear-structure effects impact negatively on the extraction of 
strange-quark matrix elements 3, 12, 13- 



Neutrino-induced reactions provide a viable alternative to parity-violating electron scat- 
tering. While the latter is primarily sensitive to the strange electric and magnetic form 
factors of the nucleon, the former is particularly sensitive to the axial-vector form factor of 
the proton — through the combination (Am — Ad — As + Ac — Ab + At). This sensitivity 
is the result of the small weak-vector charge of the proton (1 — 4sin 2 #w — 0.08) and the 
suppression of the weak anomalous magnetic moment at small Q 2 . In the above expression 
the heavy quark flavors (c, b, and t) can be eliminated using a well-defined renormalization 
group procedure ljl, 2(| . Further, the isovector combination (uj^u—dj^d) is constrained 



from neutron beta decay. This leaves the (assumed isoscalar) strange-quark contribution to 
the spin of the proton As to be determined from the elastic neutrino-proton reaction. Yet 
an absolute cross-section measurement of this reaction is an experimental challenge due to 
difficulties in the determination of the absolute neutrino flux. An attractive alternative has 
been proposed by Garvey and collaborators j^, 2^] in which the extraction of As proceeds 



through a measurement of the ratio of proton-to-neutron cross sections in neutral-current 
(NC) neutrino-nucleon scattering. This ratio is defined by the following expression: 

Rip ,n) = °JfE^EL . (!) 
o[yn — > un) 

This ratio is very sensitive to the strange-quark contribution to the spin of the nucleon as 
As [or g s A = G^(Q 2 = 0)] interferes with the isovector contribution (G^ ) with one sign in 
the numerator and with the opposite sign in the denominator [see Eq. H^. Unfortunately, 
R{p/n) is difficult to measure with the desired accuracy due to experimental difficulties 
associated with neutron detection It is for this reason that FINeSSE will focus initially 
on the neutral- to charged-current ratio (NC/CC): 

R(NC/CC) = a } Up - UP \ . (2) 

This ratio is "simply" determined from counting the number of events with an outgoing 
proton and missing mass relative to those events with an outgoing proton and a muon. Note 
that the CC reaction, being purely isovector, is insensitive to As. As such, R(NC/CC) is 
about a factor of two less sensitive to As than R(p/n). 
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From a theoretical perspective extracting As from the ratio of cross sections is also 
attractive. As a large number of the scattering events at FINeSSE will be from nucleons 
bound to a Carbon nucleus, it is important to understand nuclear-structure corrections. 
This issue came to light in experiment E374 at the Brookhaven National Laboratory (BNL) 
where it was found that 80% of the events involved neutrino scattering off carbon atoms, 
while only 20% were from free protons. As nuclear-structure corrections in R(NC/CC) 
appear to be insensitive to final-state interactions between the outgoing proton and the 
residual nucleus |23|, |24|, the ratio R(NC/CC) may be accurately computed using a much 
simpler plane-wave formalism. Indeed, in Sec. |H] we will show how the cross section ratio 
R(NC/CC) in Carbon computed in a plane- wave formalism may be expressed in a form 
that closely resembles the "Feynman-trace" approach used to calculate the cross section 
from free nucleons. We note in closing that the new generation of neutrino experiments will 
require a thorough understanding of neutrino-nucleus interactions since the detectors often 
contain complex nuclei. Experimental and theoretical work related to neutrino scattering 
from light and heavier nuclear targets may be found in Refs. H H H H m H H E2 ■ 

Our paper has been organized as follows. In Sec. |H] we briefly review the formalism 
developed in Ref. [33| for the neutral-current case and point out the main modifications 
required to make it applicable to charged- current neutrino-nucleus scattering. Our main 
results — with a focus on the sensitivity of R(NC/CC) to As — are presented in Sec. IIHI 
Finally, we summarize the main points of this work in Sec. IIV1 



II. FORMALISM 

In this section the formalism for the description of charged-current neutrino-nucleus scat- 
tering is presented. As the basic outline follows closely the neutral-current formalism devel- 
oped in Ref. (33[, we present a brief review that focusses on those modifications that arise 
from a finite muon mass. 



A. Cross section in terms of the leptonic and hadronic tensors 



The lowest-order Feynman diagram for the knockout of a bound proton via the charged- 
current reaction [v + X(Z, A) — > + p + X(Z— 1, A— 1)] is shown in Fig. ^ Here the initial 
four-momentum of the (left-handed) neutrino is k while the four-momentum and helicity of 
the outgoing muon are k' and h', respectively. The reaction proceeds via the exchange of a 
virtual W + boson with four-momentum q= (uj,q). The kinematical variables defining the 
hadronic arm are the four-momentum of the target (P) and residual nucleus (P')- Finally, 
p' and s' denote the four-momentum and spin component of the ejectile proton. Energy- 
momentum conservation demands that: 



q = k- k! = p' + P' - P . 



(3) 



The dynamical information for this reaction is contained in the transition matrix element 
given by 



iM 



7i(k', ti) 



-w 



2V2 

(pW; * f (P') 



(7 M -7 M 7 5 ) 

-ig 



2V2 



cos 9c J u {q) ^0P)>V 



(4) 
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FIG. 1: Lowest-order Feynman diagram for charged-current neutrino-nucleus scattering. 



In Eq. (jlj) the initial and final nuclear states are denoted by ^i(-P) and tyf(P'), respectively. 
Furthermore, g is the weak coupling constant, 6 C is the Cabbibo angle (cos 6c = 0.974), and 
J v (q) is the weak nuclear current operator. As only low-momentum transfers {—q^^Q 2 -C 
M^r) will be considered, the following approximation is valid 



2V(<?) 



-9vv + %(lvl M W 



M 2 W 



m 2 w 



Using this expression the transition matrix element may be written as 

M = ^^s9 c [/Z(k',/0 7 M(l-7 5 )Kk)] [{p',s';* f (P')\Hq)\MP)) 
where the Fermi coupling constant Gp has been introduced via 



G p 

71 



(5) 



(6) 



(7) 



In Eq. © /-t(k', h') is the Dirac spinor for the outgoing muon expressed in the helicity 
representation. That is, (suppressing "prime" indices for clarity), 



/i(k, h) 



Ek + m 
2E k 



( 



h- 



k 



\ 



<Mk) 



Vk 2 + m? 



(8) 



\ E k + m 

where 0h=±i(k) are two-component Pauli spinors given by 

cos(#/2) 
sin(0/2) e*+ 



>h=+l 



(k) 



i(k) 



sin(0/2) e"^ N 
cos(fl/2) 



(9) 



Here m denotes the muon mass, and 9 and are the polar and azimuthal angles of the muon 
momentum. The neutrino spinor z/(k) is directly obtained from the above expressions by 
setting the fermion mass to m = and the helicity to h —— 1. Note that left/right projection 
operators on plus/minus helicity states do not vanish in general due to the finite muon mass. 
That is, 

V L/R fi{k,h = ±l) = 1(1 T7 5 )Mk,/i = ±l) = 0{m/k) ^ . (10) 
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Further, we have adopted a non-covariant normalization for the Dirac spinors of Eq. (|5|). 

//t(k, /i)Mk, /i') = /i(k, /i) 7 V(k, tf) = <W , (11) 

a choice that is in accordance with the standard normalization of the bound-state spinor [34| 
and that is given by 

J Uiiv)U a {v)dh = l. (12) 
Following Ref. [33| the differential cross section can now be written as 



da 



Gl cos 2 6, 



° d 3 k'd 3 p'5(E k + M A - E k , - E v , - Ep^i^W^ 



2(2tt) 5 

where the leptonic tensor is given by 



-fit; 



Tr 



(13) 



(14) 



and a discussion of the hadronic tensor is postponed until the next section. 

We conclude this section with the evaluation of the leptonic tensor. To do so, both 
matrices vV and /i/Z in Eq. (fT^j) are first expressed in terms of Dirac matrices. For the case 
of the massive muon we obtain 



/u(k', h')fi(k', h'^ 



(ft 1 + m) 
2Eui 



where the four-component spin vector is given by 

1 



m 



k',E h ,\s!\ . 



(15) 



(16) 



The corresponding expression for the massless left-handed neutrino may be obtained from 
the above equations by setting the helicity to h = —1 and by taking the massless (m — + 0) 
limit. Note that in the massless limit ms M — > k ,fl . Thus we obtain, 



Kk)F(k) 



2 Eu 



5(1+7- 



(17) 



Finally, by substituting the above expressions into the leptonic tensor of Eq. (|14j). which in 
turn we separate into (// +-+ v) symmetric and antisymmetric parts, 



Ijiv — guv _|_ q\w 



A 1 



we obtain 



?r = — 
s kE J 



- (k»K' v + K'»h u - g^k ■ K ,y j 



9V» 
A 



(18) 

(19a) 
(19b) 



Note that in the above expressions the following four-vector has been introduced: 



K' = -(k' - h'ms) — ► k'S h > _i 

2 m=0 ' 



(20) 
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where the last expression denotes the massless limit. Hence, in the m^O limit the leptonic 
tensor vanishes for positive-helicity (h' = +l) but for negative-helicity (h' = —l) goes over to 
Eq. (17) of Ref. [33]. Finally, note that the following convention was adopted [35]: 

Tr (-f-f-ff-f) = > (^° 123 = -1, ^0123 = +1) • (21) 

We close this section with a comment on the conservation (or rather lack-thereof) of the 
leptonic tensor. While the antisymmetric component satisfies 

= tTl* = , (22) 

due to the antisymmetric property of the Levi-Civita tensor, this is no longer true for the 
symmetric part due to the finite muon mass. That is, 

q/^ j£ and t^q v + . (23) 



B. Differential cross section in terms of nuclear structure functions 



In Eq. ()13|) of the previous section it was shown that the differential cross section for the 
CC reaction may be written as a contraction of the leptonic tensor with the hadronic tensor, 
where the latter is defined in terms of the expectation value of the weak nuclear operator 
[see Eq. flHJ)]. That is, 



{p',s';* f (P') >(g) %(P)) (p',s'^ f (P') r(q) %(P)) 



(24) 



Although the general form of the hadronic tensor was introduced and discussed in detail in 
Ref. [33], some of its most salient features are underscored here for completeness. For the 
case of unpolarized proton emission, the hadronic tensor may be written in terms of thirteen 
independent structure functions, 

= W x g^ + W 2 q»q v + W 3 P»P U + W 4 p ,fl p' u 

+ W 5 (q^P v +P fl q 1/ ) + W 6 (q fl p >u + p^q v ) + W 7 {P fl p w + p'^P u ) , (25a) 

= ^(g^-PY) + W 9 (q»p' u - p'»q v ) + W w {Py v -p' il P v ) 

+ W u e^q a Pp + W 12 e^ aP q a p'p + W lz e^ aP P a p' p . (25b) 

Note that all structure functions are functions of the four Lorentz- invariant quantities, q^q^ = 
— Q 2 , q ■ P, q ■ p', and P ■ p'. Details on the contraction between the leptonic and hadronic 
tensors 



"puV 



l'i' " ■- ' in<" S T , (26) 

have been reserved to Appendix [XJ Yet we note that the charged-current reaction is now 
sensitive to the three structure functions W?, W*, and W%. This is in contrast to the 
neutral-current reaction (see Eq. (22a) of Ref. |33j|); the origin of this difference is the non- 
conservation of the symmetric part of the leptonic tensor due to the finite muon mass [see 
Eq. (J23)) ]. However, as the antisymmetric part of the leptonic tensor is manifestly conserved, 
both NC and CC processes are insensitive to the W 8 and Wg structure functions. 

This concludes the mo del- independent description of charged-current neutrino-nucleus 
scattering. In summary, the cross section may be parametrized in terms of eleven nuclear- 
structure functions. In principle, they could be determined by a "super" Rosenbluth separa- 
tion. In practice, however, this is not possible so we resort to a relativistic mean-field model 
to obtain explicit expressions for these quantities. This will be done in the next section. 
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C. Model-dependent evaluation of the cross section 



In the previous section a model-independent formalism was presented for charged-current 
neutrino-nucleus scattering. Specifically, the cross section was written in terms of a set of 
nuclear structure functions that parametrize our "ignorance" about the strong-interactions 
physics at the hadronic vertex. However, to proceed any further a number of approximations 
need to be made in order to obtain a numerically tractable problem. 

The first "no- recoil" approximation, detailed in Eqs. (23)-(26) of Ref. 33[, is purely 



(27) 



kinematical and yields the following expression for the angle-integrated differential cross 
section: 

= G lTy k'E k ,p'E p , sin a da I\m6d9 t W^W" 

Here a is the polar angle defining the direction of the outgoing proton having momentum 
p' = |p'| and energy E p i = \/p' 2 + M 2 . Similarly, 6 and <\> define the polar and azimuthal 
angles of the outgoing muon with momentum k' = |kH and energy = \/k' 2 + m 2 . (For 
further details we refer the reader to Fig. 2 of Ref. [33]). Finally, to compare the present 
charged-current calculation to the neutral-current one, we have integrated over the kinematic 
variables of the outgoing lepton. 

The second approximation concerns the evaluation of the nuclear matrix element 

J^tf^'-VfiP'^iq^iP)) • (28) 

First, two- and many-body components of the current operator are neglected by assuming 
that the H^-boson only couples to a single bound neutron. Second, two- and many-body 
rescattering processes are neglected by assuming that the detected proton is associated 
with the specific bound neutron to which the H^-boson had coupled to. Further, as we 
are confident that distortion effects largely factor out from the ratio of cross sections [23l 
24j . final-state interactions between the outgoing proton and the residual nucleus will be 



neglected. Finally, the impulse approximation is invoked by assuming that the weak charged- 
current operator for a nucleon in the nuclear medium retains its free-space form. That is, 

J, = tf C ~ -J%° = FriQ 2 )!, + l F 2 (Q 2 )a^ - G A (Q 2 ) W . (29) 

Here M is the nucleon mass and F\, F 2 , and Ga are Dirac, Pauli, and axial- vector nucleon 
form factors. Note that the pseudoscalar form factor has been neglected, since its contribu- 
tion is suppressed by the small lepton mass A detailed discussion of the weak charge 
current [Eq. (j2H|l ] has been reserved to Appendix El As we have assumed that the ratio of 
cross sections given in Eqs. (Q) and (J2J) are insensitive to final state interactions between 
the outgoing proton and the residual nucleus, both initial (bound) and final (free) nucleon 
propagators may be written in terms of Dirac gamma matrices — rendering the hadronic 
tensor analytical. Explicit expressions for both propagators and for the analytic (albeit 
model-dependent) hadron tensor are given in Eqs. (34-38) of Ref. |33| . 



III. RESULTS 

In this section results are presented based on the formalism outlined in Sec. |H] for 
charged- changing neutrino scattering from 12 C. The angle-integrated differential cross sec- 
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tion [Eq. ([27)1] is shown in Fig.|2]as a function of the kinetic energy T p > of the outgoing proton 
in the laboratory frame for three incident neutrino energies, namely, Ek = k = 200, 500 and 
1000 MeV. We display separately the contribution to the cross section from the lp 3//2 (solid 
and long-dashed-short-dashed lines) and ls 1//2 (dashed and dotted lines) orbitals computed 
in a relativistic mean-field approximation using the NL3 parameter set [37|. Note that be- 
cause of the finite muon mass, both negative (h' = —l) and positive (h' = +l) helicity muons 
contribute to the cross section; the two smallest contributions correspond to the positive- 
helicity case. As the energy of the incident neutrino increases, and consequently also that of 
the outgoing muon, the positive-helicity contribution (which scales as m/Ey) becomes less 
and less important until it ultimately disappears at large-enough energy. This can already 
be observed at £ fc = 500 and 1000 MeV. 




10 20 30 40 50 60 70 
T„ ( MeV ) 

16 L j 16 L 




T„ (MeV) T p , (MeV) 

FIG. 2: Differential cross section da/dE p > [Eq. (|27j) ] as a function of the outgoing proton laboratory 
kinetic energy T p i. The solid and dashed lines denote the contributions from the lp 3 ^ 2 and Is 1 / 2 
neutron orbitals in 12 C to the negative-helicity (h = — 1) cross section. The long-dashed-short- 
dashed and dotted lines are the corresponding contributions for positive helicity (h = +1). The 
three considered incident neutrino energies are £^ = 200, 500, and 1000 MeV. 

For the elementary process v + n — > + p, the threshold laboratory energy of the 
incident neutrino is approximately 112 MeV. An additional kinematical constraint that is 
strongly affected by binding energy corrections follows from energy conservation. Using the 
fact that Ey > m we obtain 

E k + (M + E B ) > m + (T p , + M) T p , < E k - E B - m , (30) 

where Eg is the (positive) binding energy of the neutron. For 12 C, the NL3 parameter set 
predicts ^(ls 1 / 2 ) ^53 MeV and E B (lp 3 / 2 ) w 19 MeV. For the particular case of a neutrino 
incident energy of E k = 200 MeV, the cross section displays a sharp cut-off for the knockout 
of the Is 1 / 2 neutron at an energy of T p /^A0 MeV. For higher incident neutrino energies, the 
maximum allowed value for the kinetic energy of the outgoing proton is already sufficiently 
large to allow the cross section to fall off smoothly to (almost) zero. Our subsequent results 
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will only focus on incident neutrino energies of 500 and 1000 MeV, as the ratio Rnc/cc 
defined in Eq. (J2J) will be measured by the FINeSSE collaboration with neutrinos in that 
energy range 0. 




o ' 

60 120 180 240 300 

^ 25 i i 
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FIG. 3: Differential cross section da/dE p > [Eq. (|27|)] as a function of the outgoing proton laboratory 
kinetic energy T p >. The target nucleus is 12 C and the incident neutrino energy is taken to be 500 
and 1000 MeV. The solid line represents the cross section calculation summed over both muon 
helicities (h = ±1) and over both bound-state (Is 1 / 2 and lp 3 / 2 ) orbitals. The dashed and long- 
dashed-short-dashed lines represent the individual contributions to the (unpolarized muon) cross 
section from the lp 3 / 2 and Is 1 / 2 orbitals, respectively. 

The cross section that results from adding the contributions from both muon helicities 
(h = ±l) and both neutron orbitals (Is 1 / 2 and lp 3//2 ) is depicted in Fig. El by the solid line. 
The dashed and long-dashed-short-dashed lines represent the calculation where we have 
summed over the two helicity values for the individual lp 3 ^ 2 and Is 1 / 2 orbitals, respectively. 
The result for the full cross section (solid line) may be compared to Fig. 8 of Ref. (38[. In 
the kinematical region in which they can be compared, there is good agreement in both the 
shape and magnitude of the cross sections. 

Next we investigate in Fig. 0] the contribution from the single-nucleon form factors to 
the differential cross section for incident neutrino energies of Ek = 500 and 1000 MeV. As 
in the previous figures, the full result is displayed by the solid line. Next in importance is 
the long-dashed-short-dashed line obtained by setting the weak Pauli form factor to zero 
(F 2 = 0). The last three lines are obtained from calculations using a single non-zero form 
factor. That is, the dashed line is obtained from the full calculation by setting Fi = F 2 = 0, 
the dotted line by setting G a = F 2 = 0, and the dashed-dotted line by setting Ga — Fi — 0. 
This figure clearly illustrates the relatively minor role played by the kinematically suppressed 
weak Pauli form factor F%. Indeed, by itself, it yields a partial cross section that both in 
magnitude and in shape shows little resemblance to the full cross section. Clearly, the two 
dominant form factors are the weak Dirac and the axial- vector form factors, with the latter 
assuming the dominant role. Yet by itself, no single form factor reproduces the full cross 
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section indicating that all interference terms, F1F2, F\Ga, and F2G a are important for the 
charged- current process. Contrast this to the neutral-current process where the Dirac form 
factor is strongly suppressed by the weak mixing angle (1— 4 sin 2 #vk~ 0.076). 




160 320 480 640 800 

T p , (MeV) 



FIG. 4: Effect of the single-nucleon form factors on the differential cross section [Eq. IJ27JI ] as a 
function of the laboratory kinetic energy T p > of the outgoing proton. The calculations include a 
sum over the two (Is 1 / 2 and lj> 3 / 2 ) neutron orbitals in 12 C and are for incident neutrino energies 
of 500 and 1000 MeV. Explanation for the various lines is given in the text. 



As mentioned earlier, systematic errors with the neutron detection make the ratio of 
neutral- to charged-current reactions R(NC/CC) a more viable alternative than the proton- 
to-neutron neutral-current ratio R(p/n). Thus, we now compare in Fig. El the cross section 
for the charged-current reaction with that for the neutral-current process: u+X(Z,A) — > 
v+p+X(Z— 1, A— 1). The solid line represents the full charged-current cross section from 
12 C, where we have summed over both bound-state orbitals and both muon helicities. The 
dashed line represents the corresponding neutral-current cross section with no strange-quark 
contribution to the spin of the proton (i.e., g s A = 0); the long-dashed-short-dashed line is 
the same calculation with g s A = —0.19. Results are shown for incident neutrino energies of 
500 MeV (top graph) and 1000 MeV (bottom graph). A comparison to Fig. 8 of Ref. [38| 
shows good agreement in both the shape and magnitude of the cross sections. The axial- 
vector form factor plays a dominant role in the neutral-current neutrino-proton reaction and 
makes this reaction particularly sensitive to the strange-quark contribution to the spin of 
the proton. Recall that the axial-vector form factor of a proton in the neutral-current case 
is given by [IF 



G A (Q 2 



9 a 9 a I 



;i-26-^) 



1.45 



9 a- 



-0.19 



(31) 



Here Gp(Q 2 ) is the axial- vector form factor of the nucleon (see Appendix El) and a value of 
g S A = — 0-19 is assumed for the strange-quark contribution to the spin of the nucleon j^]; this 
value seems to improve the agreement with the Brookhaven National Laboratory experiment 
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FIG. 5: Differential cross section da/dE p i [Eq. 1)27(1 ] for neutrinos on 12 C as a function of the 
laboratory kinetic energy T p i of the outgoing proton. The solid line represents the full charged- 
current neutrino-nucleus cross section, while the dashed (long-dashed-short-dashed) line represents 
the corresponding neutral-current neutrino-nucleus cross section using g s A = (g A = —0.19). Results 
are shown for incident neutrino energies of 500 and 1000 MeV. Note that for clarity the neutral- 
current cross sections have been multiplied by a factor of 4. 



E734 jijj . Note that this negative value of g s A leads to an increase in the proton Ga by about 
15%. 

We are now in a position to display results for the main observable of this work: the ratio 
of neutral- to charged-current neutrino-nucleus scattering cross sections R(NC/CC) defined 
in Eq. (J2J). For notational simplicity let us denote the differential cross section da/dE by a, 
where it is implied that we have summed over the Is 1 / 2 and lp 3 / 2 orbitals of 12 C as well as 
over the two values of the helicity (when appropriate). Then, because of the dominance of 
the axial- vector form factor we may write the NC cross sections as 



VNc(g s A = °) V 9 A 

As the strange-quark contribution to the spin of the nucleon is assumed to be isoscalar, the 
charged-current reaction is insensitive to it. Thus, 



R(NC/CC;g A = 0) \ g A J f A — o.w 

That is, assuming the dominance of the axial- vector form factor in neutral-current neutrino- 
proton scattering, a ~ 30% enhancement (for g% = — 0.19) in R(NC/CC) is expected from 
a non-zero strange-quark contribution to the spin of the nucleon. The R(NC/CC) ratio is 
plotted in Fig. Ifjas a function of the laboratory kinetic energy T p > of the outgoing proton for 
incident neutrino energies of 500 and 1000 MeV. The solid and long-dashed-short-dashed 
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FIG. 6: Ratio of neutral- to charged-current neutrino-nucleus ( 12 C) cross sections as a function of 
the laboratory kinetic energy of the outgoing proton T p > . The solid and long-dashed-short-dashed 
lines correspond to g s A = — 0.19 and g A = 0, respectively. The dashed line is obtained by multiplying 
the g\ = result by a constant factor of 1.32 [see Eq. (|33j) ]. The incident neutrino energies are 500 
and 1000 MeV. 



lines correspond to non-zero and zero values of g A , respectively. Further, the dashed line 
has been obtained by multiplying the g s A = result by a constant enhancement factor of 
1.32. The agreement between the solid and dashed lines indicates that the simple estimate 
given in Eq. (}3*3*j) is quantitatively correct — especially at small T p > (or equivalently small 
momentum transfer q) where the contribution from the interference term i^G^ remains 
small. While significant, the sensitivity to g A in Rnc/cc is about a factor of two less than 
in R(p/n) where both proton and neutron NC cross sections are sensitive to g A . We trust 
that after working out the systematic uncertainties in the neutron detection, this crucial 
experiment will also be performed. 

We close this section with a brief comment on the mild oscillations displayed by the 
neutral- to charged-current ratio Rnc/cc, especially at 1000 MeV. Note that this structure 
is not unique to the CC cross sections (see Fig. |HJ) but has already been observed in the NC 
cross sections of Ref. [33| (see Fig. 9). As neither the momentum distribution of the bound 
nucleons nor the nucleon form factors display any such structure, we attribute this behavior 
to a kinematical effect ori gina lly pointed out in Ref. [39| (see Fig. 1) and later reproduced 
by us in Figs. 5-6 of Ref. [33|. The mild oscillations in the single-differential cross section 
da/dT p i is a residual effect associated with the existence of a "double-humped" structure 
in the double-differential cross section d 2 a/dT p /d(cosa) (here a is the polar angle of the 
outgoing proton with kinetic energy T p i). In turn, the emergence of the double- humped 
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structure is a purely kinematical effect that results from the inability of the reaction to 
produce medium-energy nucleons. That is, high-energy neutrinos are able to produce low- 
or high-energy nucleons but not medium energy ones. While the double-humped structure is 
a robust kinematical effect, the integration over a introduces model dependences that smooth 
out — to a greater or lesser degree — some of the structure displayed by d 2 a / dT p /d(cos a) . 



IV. SUMMARY 

The distribution of mass, charge, and spin in the proton are among the most fundamental 
properties in hadronic structure. In this context, a topic that has received tremendous 
attention for over fifteen years is the contribution of strange quarks to the structure of the 
proton. In this work we have focused on the strange-quark contribution to the spin of the 
proton (g%). Elastic neutrino-proton scattering at low momentum transfer is particularly 
well suited for this study, as the axial-vector form factor of the proton — the observable 
that encompasses the spin structure of the proton [see Eq. (}3*Tf ] - dominates this reaction. 
Indeed, the two "competing" Dirac and Pauli form factors are strongly suppressed, the 
former by the weak mixing angle (1 — 4 sin 2 9w ~ 0.0 76) and the latter by the nucleon mass 
(|q|/2M). Yet in an effort to reduce systematic uncertainties related to the neutrino flux, a 
"ratio-method" has been proposed to extract g s A . Two ratios are particularly useful in this 
regard: a) proton-to-neutron yields in elastic neutrino scattering [Eq. (JTJ] and b) neutral- 
current-to-charged-current yields [Eq. (J2J)]. While the former shows a larger sensitivity to g A , 
the latter is insensitive to systematic errors associated with neutron detection. As neutrino 
experiments involve extremely low count rates, these reactions use targets that consists of a 
combination of free protons and nucleons bound into nuclei. Thus, nuclear-structure effects 
must be considered. 

In the present work we have extended the formalism developed in Ref. [33] for neutral- 



current neutrino-nucleus scattering to the charged-current reaction. In particular, cross- 
section ratios have been computed within a relativistic plane wave impulse approximation. 
Benefiting from work done by others ji^, 24]|, we justify the omission of final-state inter- 



actions by the suggestion that while distortion effects change the overall magnitude of the 
cross section, they do so without a substantial redistribution of strength. Nuclear-structure 
effects — which enter in our formalism exclusively in terms of the momentum distribution of 
the bound nucleons computed at the mean-field level — were incorporated via the accurately 
calibrated relativistic NL3 parameter set |37|. The validity of the plane- wave approxima- 
tion yields theoretical cross sections that may be displayed in closed, semi-analytic form. 
Although the structure of the weak hadronic current is the same for the neutral- and charged- 
current reactions, a few differences emerge. First, the finite muon mass results in muons 
produced with both negative and positive helicity. Further, a finite muon mass produces 
cross sections that display a sharp cut-off for low values of the incident neutrino energy. How- 
ever, for the range of neutrino energies of interest to the FINeSSE collaboration (500-1000 
MeV) the positive-helicity contribution becomes negligible. Further, while the same three 
nucleon form factors enter the neutral- and charged-current reactions, their quantitative im- 
pact differs considerably. For example, while the Dirac form factor for the neutral-current 
process is strongly suppressed (Fi (Q 2 = 0) = 0.076) it is large for the charged-current process 
(F\(Q 2 = 0) = 1). Hence, no single form factor dominates the charged-changing reaction. 
More importantly, as the strange-quark content of the nucleon is assumed to be isoscalar, 
the purely isovector CC reaction is insensitive to the strange-quark content of the nucleon. 
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This renders the ratio R(NC/CC) less sensitive to strange quark effects (by about a factor 
of 2) than the neutral-current ratio R(p/n). Still, for the value of g A = —0.19 adopted in 
this work j39|, a 30% enhancemenent in R(NC/CC) is obtained relative to a calculation 
with g A = 0. We note that our results for the charged- current cross sections were compared 
to similar calculations done in Ref. |38[ and good agreement was found in both the shape 
and the magnitude of the cross section. 

In summary, the sensitivity of the ratio of neutral- to charged-current cross sections to 
the strange-quark contribution to the spin of the nucleon g s A was investigated in a relativis- 
tic plane wave impulse approximation. The enormous advantage of this formalism is that 
our theoretical results may be displayed in closed, semi-analytic form. The central moti- 
vation behind this work is the proposed FINeSSE program that aims to measure g s A with 
unprecedented accuracy via the neutral- to charged-current ratio R(NC/CC). By adopting 
a value of g s A = —0.19, an increase in this ratio of approximately 30% was found relative to 
the g A = result. While sensitive, this is less so than the corresponding ratio of proton- 
to-neutron yields R(p/n) in neutral current neutrino-induced reactions. This measurement, 
however, has been hindered by difficulties associated with neutron detection. We trust that 
this difficulty may be overcome so that this crucial program may get off the ground. 
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APPENDIX A: LEPTONIC-HADRONIC CONTRACTION 



In Sec. Ill Al it was shown that the charged-current cross section could be expressed as the 
contraction of a leptonic tensor [Eq. (|14j)] with a hadronic tensor W^ v [Eq. (j^ljl ] written 
in a model-independent way in terms of thirteen independent structure functions. In this 
appendix we carry out the contraction, which we separate into symmetric and antisymmetric 
parts. That is, 

^W, v = l%W^ + l%W^ . (Al) 
Here the symmetric part is given by 

(db) ^ vW ° v = {- W ^ k - K ') + w * fib) + w * M p ) + w * a(po 

+ W 5 f 2 (P,q) + W 6 f 2 (q,p') + W 7 f 2 (P, P ')') , (A2) 

while the antisymmetric part by 



1 



C W a = * [Wio eT a %K' v P^ + W n fc(q, P) + W 12 f 3 (q, p') + W 13 fc(P, p') 

(A3) 



Note that the following four-vector has been defined: 

K' = -(k'- h'ms) — ► k'5 h i _i . (A4) 

2 m=0 ' 

Further, for simplicity the following three functions have been introduced: 

fx(x) = 2(k ■ x){K' ■ x) - x 2 (k ■ K') , (A5a) 
f 2 {x, y) = {k- x){K' -y) + (k- y){K' ■ x) — (x ■ y)(k ■ K') , (A5b) 
h{x,y) = (k ■ y){K' ■ x) - (k ■ x)(K' ■ y) . (A5c) 

From Eq. ()A2|) we see that the three structure functions W 2 , W 5 and W 6 do contribute 
to charged-current neutrino-nucleus scattering, in contrast to the neutral-current case (see 
Eq. (22a) of Ref. [33j]). This is due to the lack of conservation of the symmetric part of the 
leptonic tensor as a result of the finite muon mass [see Eq. (jl9aj) ]. Note, however, that in 
the massless limit fi(q) = f 2 (P,q) = f2(q,p') = 0, as required. Finally, due to the form of 
Eq. f!19b|) . the charged-current process remains insensitive to the two structure functions W& 
and Wg. 



APPENDIX B: SINGLE NUCLEON FORM FACTORS 

In Sec. Ill CI it was shown that in the impulse approximation, the single nucleon current 
probed in the charge- changing reaction may be written in the following standard form: 

J, = J»° ~ J™ = F.iQ 2 )!, + iF2{Q 2 )°^ - G A (Q 2 h, l5 , (Bl) 

where Fi, F 2 , and Ga are the Dirac, Pauli, and axial-vector form factors, respectively and the 
pseudoscalar form factor has been neglected. To understand the structure of the vector form 
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factors (Fi and F%) we invoke the conservation of the vector current (CVC) hypothesis. To 
start, one parametrizes the nucleon matrix elements of the isovector electromagnetic current 
in the following standard form: 

(N(p',s',t')\j™(T=l)\N(p,s,t)) = (N(p',s',t')\qj^q\N(p,s,t)) 

U(p', s>) [if + zF«(Q>^] U(p, s) (r 3 ) ^ , (B2) 

where q = (u, d) is an isospin doublet of quark fields, and F± and are the isovector 
Dirac and Pauli form factors of the nucleon, respectively. In turn, these are given in terms 
of proton and neutron electromagnetic form factors as follows: 

Fi 1] {Q 2 ) = ~ (// P) (Q 2 ) - Ft\Q 2 )) , (i = 1,2) . (B3) 

The CVC hypothesis is a powerful relation that assumes that the vector part of the weak 
charge- changing current may be directly obtained from the isovector component of the elec- 
tromagnetic current. That is, 

JEM (T = 1} = y(3) = JCC (±) = y(l) ± • yf) = ^ q . (B4) 

Thus, a determination of the electromagnetic form factors of the nucleon — which has been 
done experimentally — fixes the vector part of the charge-changing currents to: 

(N(p',s',t')\J^(±)\N(p,s,t)) = U(p',s') [Fl l \Q 2 ) 1 » + iF?\Q 2 )a^]u{p,s) (2r±)^ . 

(B5) 

In this way the vector form factors of Eq. (jBljl are then simply given by 



F l (Q 2 )=2F^(Q 2 )=F t {p \Q 2 )-F^(Q 2 ), (< = 1,2). (B6) 



Paraphrasing Ref. |41j : CVC implies that the vector part of the single-nucleon matrix element 
of the charge- changing weak current, whatever the detailed dynamic structure of the nucleon, 
can be obtained from elastic electron scattering through the electromagnetic interaction! 

A similar procedure may be followed to determine the axial-vector form factor Ga in 
terms of the isovector axial- vector current. That is, 

^ C 5 C (±) = ± i Af = g 7 , 75 (^P) q , (B7) 



so that 



(N(p', S \t')\J^ c (±)\N(p,s,t)) = G A (Q 2 )U(p>,s>) W U(p,s) (r±) tt; . (B8) 

As before, the above expression neglects the contribution from the pseudoscalar form factor. 

We finish this section by parameterizing the various nucleon form factors in terms of their 
known Q 2 = values times form factors of a dipole form. This is identical to the procedure 
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employed in Appendix A of Ref. |33j|) for the neutral-current reaction. We obtain 

F[ P \Q 2 ) = ( 1+ ^ Ap) ) , F^(Q 2 ) = GUQ 2 ) , (B9a) 

F[ n \Q 2 ) = (^y^ 1 ) G V D (Q 2 ) , Ft\Q 2 ) = ( An( 1 1 ^ 7 " ?) ) G ^ (g2) ' (B9b) 

G A (g 2 ) = ^(g 2 ), (B9c) 

where a dipole form factor of the following form is assumed: 

G V D (Q 2 ) = (1 + Q 2 /M 2 )~ 2 = (1 + 4.97r)- 2 (BlOa) 
G£(Q 2 ) = (1 + Q 2 /M 2 )- 2 = (1 + 3.31r)- 2 (BlOb) 
r/ = (1 + 5.6 t)- 1 t = Q 2 /(AM 2 ). (BlOc) 

Finally, for reference we display the value of the various nucleon form factors at Q 2 = 

Fj p) (0) = l, F 1 (n) (0)=0, (Blla) 
F 2 (p) (0) = A p = +1.79 , F 2 (n) (0) = A n = -1.91 , (Bllb) 
G A (0) = 9A = +1-26 . (Bile) 
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